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Abstract 

We cloak a region from a known incident wave by surrounding the region with 
three or more devices that cancel out the field in the cloaked region without 
significantly radiating waves. Since very little waves reach scatterers within the 
cloaked region, the scattered field is small and the scatterers are for all practical 
purposes undetectable. The devices are multipolar point sources that can be 
determined from Green's formula and an addition theorem for Hankel functions. 
The cloaking devices are exterior to the cloaked region. 
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1. Introduction 

Interest in cloaking has surged, as reflected in the many recent reviews [|3l[T0l 
|2). We introduced a new kind of cloaking for the two-dimensional Helmholtz 
equation lfTTl[T2) . Our approach uses active sources (cloaking devices) to hide 
objects placed in an external region. The advantages of our approach are: (a) by 
the superposition principle a cloak can be designed for a broad band of frequen- 
cies (excluding discretely many frequencies where the object being cloaked, if 
non-absorbing, "resonates") and (b) the cloak does not need materials with ex- 
treme properties which are hard to realize and dispersive, as is the case in most 
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transformation based cloaking strategies (see e.g. (H |2TJ [T5l [6l OH - though an 
exception is U6J). A significant drawback of our approach is that we assume 
full knowledge of the incident field. Also the active sources contain a monopole 
term which may be problematic for applications in electromagnetism. 

The problem of finding source distributions for cloaking is clearly ill-posed 
in the sense that if it admits one solution then it admits infinitely many solutions. 
In IfTTl [T2ll we computed particular solutions involving three point- like devices 
by solving a constrained least-squares problem with the singular value decom- 
position (SVD). In an effort to explain rigorously our previous results, we use 
the Green representation theorem for the Helmholtz equation (in short Green 's 
formula, see e.g @), to derive explicitly a particular solution in terms of the in- 
cident field (Theorem[TJ). Another cloaking method based on Green's formula is 
the active interior cloak introduced by Miller iTTTll . which uses single and dou- 
ble layer potentials to cancel the incident field inside a closed curve. The same 
principle is used in active sound control E21 I71I to suppress noise. Here we apply 
an addition theorem for Hankel functions to replace the source distribution on 
a curve by a few active sources, effectively connecting the cloaked region with 
the exterior. The sources we get are, as the ansatz used in IfTTl fT2ll . multipolar 
point-like sources: they are given as a series of cylindrical radiating solutions to 
the Helmholtz equation, centered at a few points. We also establish convergence 
of the series to the fields required for cloaking and give a specific configuration 
of sources (§|2.3[) similar to the one we found empirically in IfTTl [T2ll . 



Other methods for obtaining exterior cloaks include those based on comple- 
mentary media lfT3ll . surface plasmonic resonances Il23ll . anomalous resonances 
in the vicinity of a superlens ffTj^ |20l [19) and waveguides Il24l . 

A different idea is that of illusion optics [fT4| where the goal is to hide an 
object and make it appear as another object. The Green's formula based approach 
that we present here also explains why this can be done using active devices, as 
was recently observed numerically lf25fl . We give a way of explicitly constructing 
the devices to such effect, without the need for solving a least- squares problem 
(see Remark [4]). 

We work in the frequency domain at a fixed angular frequency co. In a 
medium with constant speed of propagation c, the wave field w(x, to) satisfies 
the Helmholtz equation 

Au + k 2 u = 0, forxeM 2 , (1) 



where k = In/ A is the wavenumber and A = 2nc/u) is the wavelength. For 
simplicity we drop the dependency on the frequency and write w(x) = u(x, oS). 
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Remark 1. We assume that the frequency to is not a resonant frequency of the 
scatterer we want to hide. Resonant frequencies are left for future studies. 

Remark 2. Calling a cloak "active" can be ambiguous as one can refer to a 
cloak that can hide active sources ^ or a cloak that uses active sources to hide 
objects KT7\ 1771 [721/ . The cloaking method we present here is "active " in both 
senses as we use active devices to hide objects and in some simple situations 
it is possible to hide sources ( see Remark \3ty. However here we focus only on 
cloaking scatterers. 



2. Green's formula cloaks 

We present the active interior cloak [[T7l and using an addition theorem for 
Hankel functions show how this cloak can be replaced by a few multipolar 



sources (§2.2). The price to pay for using a finite number of sources is that 
the cloaked region and the region from which the object is invisible are smaller 
compared to the active interior cloak. We then make some geometric consider- 



ations (§ |2.3[ ) to show that with the particular approach of § |2.2| we need three or 
more sources to get a non-empty cloaked region. 

2.1. Active interior cloak 

Denote by D the region of IR 2 that we wish to cloak from a known incident 
field Uj. We also assume from now on that D is a simply connected bounded 
region of class C 2 . The arguments in ^2] can be easily generalized to the case 
where D is composed of several simply connected components. In order to cloak 
D, we construct a solution u c i to the Helmholtz equation ([T]) (in IR 2 excluding the 
boundary of D) such that 

(-u;(x) for x 6 D 
,V (2) 
otherwise. 

Hence the total field uj + Ui vanishes in D and is indistinguishable from U[ out- 
side D. If a scatterer is placed inside D the scattered field u s resulting from the 
incident field u c i + u t is zero. Assuming u t is an analytic solution to the Helmholtz 
equation ([I]) inside D, a field uj satisfying can be constructed using Green's 
formula (see e.g. [5]) 

u d (x)= [ J5 y (-(n(y)-V yWi (y))G(x,y) + a! (y)n(y)-V y G(x,y)), (3) 

JdD 
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Figure 1: The field uj generated by Green's formula (|3]l The region D is the disk in thick lines 
with radius 5A and the incident field is a plane wave with angle 57r/13. The field is very 
close to the incident field inside D and zero outside D. The integral in Q is evaluated with the 
trapezoidal rule on 2 8 equally spaced points on dD. The axis units are in wavelengths A. 

where n(y) is the unit outward normal to D at a point y on the boundary dD and 
the Green's function for the two dimensional wave equation is 

G(x,y) = ^\k\x-y\), (4) 

where H„ is the n-th Hankel function of the first kind [1, §9]. The first term 
in the integrand in ([3]) can be interpreted as the potential due to a distribution of 
monopoles on the boundary dD (the single layer) while the second term can be 
interpreted as the potential due to a distribution of dipoles oriented normal to dD 
(the double layer). 

In the frequency domain, the cloaking scheme we obtain is the same proposed 
by Miller [fTTTl . where the single and double layer potentials in Q are simulated 
with many sources completely surrounding the cloaked region D. We give an 
example of the field u c i generated by Green's formula ([3]) in Figure [T] The effect 
of this active interior cloak can be seen in Figure [2j where it is used to hide a kite 
shaped scatterer 01 with homogeneous Dirichlet boundary conditions. The field 
is virtually zero inside the cloaked region, so there is very little scattered waves 
to detect the object. 

2.2. Active exterior cloak 

Our aim is to replace the single and double layer potentials on dD appearing 

in ([3]) by a few ridev multipolar sources (what we call "cloaking devices") located 
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Figure 2: Active interior cloak (a) inactive and (b) active. The region D where Green's formula is 
applied is the circle of radius 5 A in thick lines. The field inside the cloaked region D is virtually 
zero, so that outside the cloaked region the field is indistinguishable from the incident plane wave 
with direction 5n/13. The axis units are in wavelengths A. 

at some points x 7 $ dD. The advantage being that the cloaked region is no longer 
completely enclosed by a surface. The field generated by such sources can be 
written formally as 

"dev °° 

"i ex l) (x) = Yj E b jM x - X A ( 5 ) 

7=1 m=-oo 

where the coefficients bj <m 6 C are to be determined and 

y m (x) = //( i 1 ^|x|)exp[/marg(x)] 

are radiating cylindrical waves. Here arg(x) denotes the counterclockwise ori- 
ented angle from the vector (1,0) to the vector x. By a radiating solution to 
the Helmholtz equation, we mean it satisfies the Sommerfeld radiation condition 
(see e.g. El). 

In IfTTl [T2ll we presented a numerical scheme based on the singular value 
decomposition (SVD) to compute the coefficients bj %m in a way that w^ ext) approx- 
imates Ud as defined in ([2]). We give next in Theorem [T] one way of obtaining 
these coefficients explicitly such that u\ (%) = aj(x) for x in a certain region 
R c IR 2 , effectively reducing the cloaked region to D n R. 

Let us assign to each source x ; a segment dDj of the boundary. These seg- 
ments are chosen such that they partition dD and dDj n dDj is empty or a single 
point when i £ j. An example of this setup is given in Figure [3j The coefficients 
bj m are given next. 
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Theorem 1. Multipolar sources located at Xj $ dD, j = 1, . . . , n c i ev , can be used 
to reproduce the active interior cloak in the region 



ndev 



xei 2 



|X - Xy| > SUp y - Xj 



7=1 I yedD J ) 

The coefficients bp n in ([5]) such that u^ xt \x) = Ud(x)for x e R are given by 

bj, m = I dS y {(-n(y) • V y « ; (y)) U m (y - x •) 

Jaoj v ^ (6) 

+ K/(y)n(y) • V y C/ m (y - x 7 )} 

/or j = 1, . . . , n ext and for m e Z. Here U m (x) are entire cylindrical waves, 

U m (x) = 7 m (£|x|)exp[z'marg(x)]. 

Proof. Clearly, the integral over the whole boundary 3D in Green's formula Q 
can be written as the sum of the integrals over the segments dDj. For segment 
dDj, we use Graf's addition formula 01 §9.1.79] to express the Green's function 
G(x, y) as a superposition of multipolar sources located at Xj, 

G(x,y) = ^H ( l \k\x-Xj-(y-Xj)\) 

i ^ 0) 

= T 2j y m(X-Xy)f/ m (y-X ; ), 

m=—co 

where the series converges absolutely and uniformly on compact subsets of |x - Xj\ > 
|y - Xj\ (this can be seen by e.g. adapting Theorem 2.10 in [4] to two dimen- 
sions). Splitting the integral in Green's formula ([5]) into integrals over dDj and 
using the expansion (|7]) we get for x e R, 



"dm 

■PW = E 



£ I dS y {(-n(y) • V y «,-(y)) ^ V m (x - x 7 )t/ m (y - x 7 ) 

j=\ ^ dD J m=-co 

(8) 

+ a ; (y)n(y) • V y ^ V m (x - Xj)U m (y - x y )}. 



The desired result (|6]) can be obtained by rearranging the infinite sum and the 
integral. For the first term in the integrals in ([8]), the uniform convergence of the 

6 



series ([?]) for y e dDj allows us to switch the infinite sum and the integral. The 
second term involves the gradient 



V y C/ m (y - xj) = k- -J' m (k |y - x ; |) exp[-zm arg(y - xj)] 



y - x ./ 



(9) 



+ J m (k |y - Xj\)(-im exp[-z'm arg(y - x ; )])V y arg(y - xj). 
Since we assumed x ; t 3D, the gradient 



„ , , (y-x 7 ) x x 

V y arg(y - x 7 ) = -, with x = 

y - x / 



-X2 



(10) 



is bounded for y e dDj. Using the series representation for Bessel functions (see 
e.g. HI §9.3.1 ] and flU §3.4]) we can get the estimates 



t n-\ 



J n (t) 



2 n n\ 



(l+(9(l/n)), 



2 n (n- 1)! 



(1+0(1//*)) 



,.„._ 2"(i»-l)! 
nit" 



(11) 



(l+0(l/n)) 



valid for t > as — » oo, uniformly on compact subsets of (0, oo). Using ( |TT| ) 
and the expression for the gradient (|9]) we can estimate the terms 

im-l \ 



V m (x - Xj)V y U m (y - xj) = O 



y - x 7 



X - X; 



+ 



y- x ; 



m- 



X - X; 



(12) 



as m — » oo uniformly on compact subsets of |x - x ; | > |y - x ; |. Therefore the 
series in the second integrand of ([8]) converges absolutely and uniformly in dDj 
and the infinite sum and the integral can be switched. □ 

Note that Theorem [T] does not guarantee that the effective cloaked region 
DC\R is not empty. However we do have that the device's field vanishes far away 
from the devices (i.e. i^ ext) (x) = for |x| large enough) because M. 2 \R is bounded. 
Later in §2.3 we give a specific configuration where D P\R ^ 0. 

Remark 3. In order to guarantee that the field Ud(x) in Green 's formula (|3]) 

vanishes outside D, we need an analytic incident field m,(x) inside D. If the field 

Uj(x) is a radiating and C 2 solution to the Helmholtz equation ([T]) outside D (as 

is the case when there are sources and non-resonant scatterers inside D), then 

Green's formula ([3]) converges outside D to Ui(x) and inside D to zero (see e.g. 

ft5}l). This is the principle behind noise suppression \7\ and could be used to 

cloak active sources and scatterers in D, assuming they are known. 
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Figure 3: Active exterior cloak construction. The contribution of portion dD y to the single and 
double layer potentials in Green's formula ([3]) is replaced by a multipolar source located at Xj t 
3D. 

Remark 4. Clearly the same Green 's formula approach can be used to create 
illusions with active sources 125]: one can simultaneously cloak an object and 
generate waves that correspond to the scattering from a completely different ob- 
ject ( the "virtual object"). All we need is knowledge of the scattered field ul"''(x) 
generated by hitting the virtual object with the incident Ui(x). Since u v s irt (x) is a 
radiating solution to the Helmholtz equation, to achieve illusion simply subtract 
u v s ' rt from Ut in ([6]) (this is assuming that u v s irt (x) is C 2 outside D). 

Remark 5. The proof of Theorem^generalizes easily to the Helmholtz equation 
in three dimensions. We leave this generalization for future studies. 

2.3. An explicit example of an active exterior cloak 

In this Section we describe one possible realization of the cloak configuration 
presented in Theorem [T] First notice that with this particular method we need to 
have ndev > 3 for a non-empty effective cloaked region D n R (with D and R 
given as in Theorem[T]). This is consistent with our numerical results in [fTTT IT2T1 . 
as we observed that at least three devices are apparently needed to cloak plane 
waves with an arbitrary direction of propagation. 

To see that at least three devices are needed, first notice that M 2 \i? = lj"=i Bj, 
where the Bj are disks centered at the j'-th device location x 7 and dDj c Bj. We 
thus get 

dD cQ Bj. 

7=1 
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If we have only one device (n^ = 1), then D <zM}\R and the effective cloaked 
region is empty. If we have two devices and D is simply connected, then we must 
also have D c M?\R, as the union of two non-disjoint disks is simply connected. 

With three devices we take as an example the configuration shown in Fig- 
ure [4] Here the devices are located at a distance 6 from the origin and are the 
vertices of an equilateral triangle. The region D where we apply Green's for- 
mula is the disk of radius cr centered at the origin. The circle 3D is partitioned 
into three arcs 3D j = 1 , 2, 3 of identical length which are chosen so that the 
distances sup y&9D |y - x ; | are equal for j = 1, 2, 3. 

Simple geometric arguments show that the region R of Theorem [j] is the 
complement of the union of the three disks in gray in Figure [4] with radius 

r((T, 6) = ((cr - 6/2) 2 + 36 2 /4) i/2 , (13) 

and centered at Xj, j = 1, 2, 3. To get an idea of the dimensions of the effective 
cloaked region R n D (in green in Figure [4]), we look at the radius of the largest 
disk that can be inscribed inside. This disk has radius 

r eS (o-,S) = 6-r(o-,S). (14) 

Thus for fixed 6 the largest effective cloaked region is obtained when cr = 6/2, 
which corresponds to the case where the intersection of two of the disks in gray 
in Figure [4] is a single point. Then the radius of the largest disk that can be 
inscribed inside R n D is 

r * eS (S) = (1 - V3/2)<5 « 0.135. (15) 



3. Numerical experiments 

We compare the Green's formula based method with the geometry described 



in 



2.3 and the maximal cloaked region size (i.e. cr = 6/2) to an SVD based 
method [fTTl[T2l . In the SVD based method three distances are needed to describe 
the devices and the cloaked region: the distance 6 from the devices to the origin, 
the radius a of the cloaked region and the radius y of the circle where we enforce 
that the device's field be small. The numerical experiments in [fTTl[T2l were done 
with a = 6/5 and y = 26. To compare the cloaked regions for both methods we 
choose 

a = r* eS (6) = (1 - Y3/2)5 
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Figure 4: A configuration for a cloak with three devices forming an equilateral triangle such 
that \xj\ = 6, j = 1, 2, 3. The region D where we apply Green's formula is the disk of radius <x 
centered at the origin. The region R 2 \R appears in gray. The green colored region is the effective 
cloaked region D C\R. 

and leave y unchanged. In this way the cloaked region for the SVD method 
should be the largest disk that fits inside the effective cloaked region D C\R of the 
Green's formula method. 

We first compare in Figure [5] the device's fields and the total field in the 
presence of a scatterer for both the SVD method and Green's formula method. 
Here the devices are at a distance 5=10/1 from the origin, with k = 1 and A = In. 
In all our numerical experiments the series Q is truncated to m = -M, . . . , M. 
For the SVD method we follow the heuristic M{8) = \(kd/2)(l + V3/2)] in JT2|, 
which for the setup of Figure [5] gives M(10A) = 59. For comparison purposes we 
used the same number of terms in the Green's formula method. In both methods, 
the device's field cancels out the incident field in region near the origin without 
changing the incident plane wave. The region where the total field vanishes is 
larger for the Green's formula method than for the SVD method. And for the 
former method, the cloaked region seems larger than what is predicted by § |2.3[ 

The fields near the devices (the "urchins" in Figure [5]) are very large as can 
be expected from the asymptotic behavior of Hankel functions at the origin 

H<P(t) = 0(t~ lnl ) as t -> for n e Z\{0}. 

The Green's formula allows us to replace each of the devices by a closed curve 
containing the device and with appropriate single and double layer potentials. 
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Figure 5: Comparison of the Green's formula method (left column) and the SVD based method 
ifTTirrZll (right column). The circle in thick lines is where we apply Green's formula. The dashed 
circle is where we enforce that the field be close to the incident field and the solid circle is where 
we enforce that the device's field is close to zero. The first row corresponds to the device's field 
and the second row to the total field. 



The curves could be chosen as circles outside of which the device's field has 
reasonable values (e.g. less than 100 times the magnitude of the incident field). 
For both methods these circles do not touch, leaving "throats" connecting the 
cloaked region to the exterior, so the cloaked region remains outside these "ex- 
tended" devices. 

We also considered larger cloaked regions in Figure |6} keeping the same 
configuration as in Figure [5] but taking 6 e [5,50]/l. To evaluate the cloak per- 
formance we plot in Figure [6] the quantity 

lk + ^ lb , on the circle |x| = (1 - V3/2)5, 

Ik-lb 
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which measures how well the device's field cancels the incident field inside the 
cloaked region and 

WUdWtf -Ill or 

- — - — , on the circle |x| = 26, 

Wi\\ L 2 

which measures how small the device's field is relative to the incident field far 
away from the devices Q With the same number of terms M{6) the SVD method 
outperforms the Green's formula method in both the cloaked region and far away 
from the devices. When 2M{6) terms are used for the Green's formula method, 
the relative errors improve, but the convergence of the series in ([5]) seems slower 
in the cloaked region than far away from the devices. 

We estimate in Figure[7Jthe size of the "extended" devices (i.e. the "urchins" 
where the device's field is large) by assuming they are disks centered at the de- 
vice's locations x,. The disk radius for the device centered at x, is estimated by 
finding the closest intersection of the level set |wd(x)| = fi with each of the seg- 
ments between x, and the other two devices. The quantity plotted in Figure [7j 
is the maximum of these distances rescaled by 6 and is always below the radius 
at which the extended cloaking devices do not leave gaps with the exterior (for 
both cut-off values f3 = 5 and 10). The cost of having more terms in Green's 
formula is that the extended devices leave narrower throats. Presumably in the 
limit M — » oo these extended devices touch and correspond to the region E?\R 
in Theorem Q] 

Remark 6. A natural question is whether changing the integrals over dDj in ^ 
to integrals over subsets of dDj gives good cloaking devices. Since in the region 
R of Theorem [7J the device 's field is identical to that of the active interior cloak, 
we can reformulate the question as follows: does making small openings in 3D 
give a good active interior cloak? This is not the case because the resulting 
active interior cloak has fields identical in R to the fields for the cloak taking 
all of 3D into account minus the fields obtained by integrating Green 's formula 
([3]) on the portions of 3D that were excluded. The excluded portions have a 
monopole and dipole distribution that radiates and spoils the cloaking effect 
that we are after, even for small openings. 



'The values of ||Mrf||/||M;|| reported in Figure 4b in fTZl are, due to a normalization mistake, up 
to a factor of VlO smaller than what they should be. In the logarithmic scale we use to display 
this quantity, the resulting shift is small and our conclusions remain the same. 
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Uj + u d \\/\\ui\\ on |x| = (1 - V3/2)5 IWI/INI on |x| = 26 



10' f 10° 




(a): 6 (in A) (b): 6 (in A) 

Figure 6: Comparison of the cloak performance for the S VD method with M(5) terms (blue) and 
the Green's identity method with M{8) terms (red) and 2M(5) terms (green). In (a) we measure 
how small is the total field inside the cloaked region and in (b) how small is the device field far 
away from the devices. 




Figure 7: Estimated device radius relative to 6 for different values of 6. The dotted line corre- 
sponds to the radius of the devices if they were touching circles. The SVD method with M(6) 
terms is in blue, the Green's identity method with M(5) terms is in red and with 2M{6) terms in 
green. 
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